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Abstract 

In  this  report  the  concept  of  a  measurement  technique  for  calculating  sound 
intensity  in  the  frequency  domain  is  discussed  and  also  how  such  a  measurement 
system  can  be  implemented  in  practice  by  using  a  frequency  domain  analyser. 

The  technique  described  employs  a  dual-channel  FFT  analyser  to  obtain  a  cross 
power  spectrum  from  the  two  microphones  from  which  socmd  intensity  is 
calculated  .  This  approach  enables  a  general-purpose  FFT  aiudyser  together  with  a 
micro-computer  to  perform  the  function  of  a  dedicated  sound  intensity  analyser. 
The  application  of  sound  intensity  measurement  technique  in  sound  power 
determination  of  a  reference  source  is  presented. 
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Measurement  of  Sound  Intensity  and  Sound 

Power 

1.  Introduction 


This  report  describes  a  new  and  practical  method  for  measuring  and  analysing 
the  noise  signature  caused  by  air>bome  and  structure-bome  noise.  If  the  soniul 
intensity  of  noise  sigttatiires  can  be  measured,  it  is  then  possiUe  to  identify  aiui 
minimize  the  average  acoustic  spectrum  and  emission  of  the  Navy’s  surface  shfys 
and  subnrarines. 

Conventional  techniques  for  measuring  sound  power  [ref3,6],  or  the 
transmission  loss  of  sound-absorbing  iruiterials,  require  special  facilities  sudt  as 
anechoic  or  reverberant  chambers.  This  report  deals  with  a  measurement 
technique  for  sound  intensity  that  is  an  attractive  alternative  to  these  conventional 
techniques. 

Measuring  sound  power  by  means  of  soimd  intensity  has  a  munber  of 
advantages  including  the  ability  to  locate  the  sources  and  sinks  of  sound,  to 
determine  sound  power  from  a  source  in  situ  even  in  a  noisy  rruchinery  room, 
and  the  ability  to  map  the  intensity  of  the  sound  source.  These  carmot  be  done  by 
the  traditional  techniques  except  under  controlled  acoustic  envinmment. 

H.F.  CMson  laid  down  the  theoretical  background  for  intensity  measurements  in 
1932  [ref.l]  but  it  was  not  until  the  1970s  that  the  electronic  instruments  required 
for  a  reliable  measurement  of  sound  intei»ity  became  available  [ref5].  A  sound 
intensity  measurement  system  comprises; 

1.  An  intensity  probe  for  sensing  the  sound  signal 

Currently  there  are  two  categories  of  probe  in  widespread  use.  The  ”p-u”  probe 
combines  a  pressure  transducer  with  a  particle  velocity  transducer.  The  ’’p-p” 
probe  is  made  up  of  two  matched  microphones  (pressure  transducers)  separated 
by  a  spacer  for  measuring  sound  pressure  and  particle  velocity  [ref.4].  This  report 
discusses  the  use  of  the  "p-p"  type  of  probe,  upon  which  the  measurement  system 
is  based. 


2.  An  analyser  for  processing  the  sound  signal 


The  analyser  may  be  one  of  two  types.  It  may  have  either  digital  or  analog  filters, 
and  operate  in  the  time  domain  (as  do  Type  4437  and  2134  Bruel  &  Kjaer  sound 
intensity  analysers);  or  it  may  be  a  dual  channel  FFT  analyser  which  operates  in 
the  frequency  domain  of  the  sound  signals  [ref.lO]. 

The  former  are  faster,  have  real-time  processing,  and  produce  output  in  the 
octave  or  fractional  octave  frequency  bands  that  are  frequently  used  in  acoristic 
measurements.  However,  if  the  analysis  function  can  be  performed  on  a  geiwral 
purpose  FFT  analyser,  it  has  the  important  advantage  of  eliminating  the  need  for 
an  expensive  and  limited  piece  of  equipment.  This  report  describes  how  a  dual 
chaimel  FFT  analyser  can  be  used  to  determine  sound  intensity,  by  taking  the 
cross  power  spectrum  from  two  microphones. 


2.  Sound  intensity  and  its  applications 


Sound  intensity  describes  the  rate  of  energy  flow  (ie.  the  power  flow)  per  unit 
area  at  a  point  in  space.  In  the  SI  unit  system,  the  unit  for  Sound  Intensity  is  watts 
per  square  metre. 

In  contrast  to  sound  pressure  which  is  a  scalar  quantity,  sound  intensity  is  a 
vector  quantity  as  it  has  both  direction  and  magnitude.  Usually,  sound  intensity  is 
measured  in  the  direction  normal  to  the  surface  through  which  there  is  a  net  flow 
of  sound  energy. 


2.1  Sound  Power  and  Sound  Pressure 

Sound  power  can  be  considered  as  the  stren^  of  the  source  whidi  gives  rise  to  a 
sound  pressure  in  space.  The  cause  of  noise  is  sound  power;  what  we  hear  is  its 
effect,  sound  pressure. 

Sound  power  is  virtually  independent  of  the  environment.  A  flve  watt  source 
will  always  give  an  energy  output  of  five  watts;  it  does  not  matter  if  the  source  is 
placed  in  a  Imge  or  small  room,  or  if  there  is  another  source  present. 

Sound  pressure,  on  the  other  hand,  is  a  quantity  that  does  depend  on  the 
acoustic  envirorunent;  on  the  distarKe  from  tlw  point  of  measurement  to  the 
source,  on  the  size  of  the  room  in  which  the  sound  source  is  placed,  on  the  sound 
absorption  coefficients  of  surrounding  surfaces,  and  so  on  [ref .6] 

Because  of  its  independence  from  the  acoustic  environment,  sound  power  is  a 
good  descriptor  of  the  source  strength  whereas  sound  pressure,  a  quantity  that 
relates  directly  to  hearing  damage  and  noise  armoyance,  is  a  good  descriptor  for 
the  effects  of  sound  on  people. 
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2.2  Applications  of  Sound  Intensity  Measurements 
2.2.1  Sound  Power  Determination 


In  determining  sound  power  radiated  from  a  source,  a  surface  enclosing  the 
source  is  defined  as  the  control  surface.  At  this  surface,  the  sound  intensity  is 
measured  in  the  direction  normal  to  the  surface  and  the  total  flow  rate  of  energy 
outward  is  determined. 

Before  the  development  of  measurement  techniques  for  sound  intensity,  the 
sound  power  of  a  source  was  usually  calculated  Aom  sound  pressure 
measurements.  These  measurements  need  to  be  made  within  special  structures 
with  known  acoustic  properties,  such  as  anechoic  or  reverberant  chambers.  Such 
structures  are  quite  costly,  and  the  sound  source  being  measured  must  first  be 
removed  from  its  normal  working  environment  [ref.5,6]. 

On  the  other  hand,  sound  intensity  can  be  measured  in  virtually  any 
environment.  This  allows  the  measurement  to  be  done  on-site.  By  using  sound 
intensity  measurements  on  a  control  surfa<%  as  described  above,  the  sound  power 
of  a  given  noise  source  can  be  determined  even  in  the  presence  of  other  radiating 
sources. 

Theoretically,  background  noise  will  have  no  effect  on  the  measurement  of 
sound  power  because  sound  intensity  is  a  vector  quantity.  Background  noise  is 
caused  by  soiuces  outside  the  control  surface,  and  so,  according  to  Gauss' 
theorem  [ref. 2],  the  net  flow  of  background  noise  through  the  control  surface  is 
zero.  (This  applies  only  to  a  stationary  background  noise  and  where  there  is  no 
sound  absorption  materials  inside  the  control  surface.) 

In  practice,  a  sound  field  generally  consists  of  active  part  (free  field  plane 
wave),  reactive  part  (standing  wave)  and  diffuse  part  (reverberant  field)  [ref.l2] . 
A  sound  intensity  measurement  system  only  responds  to  the  active  part.  When 
sound  intensity  is  measured  in  the  presence  of  a  strong  extraneous  background 
noise  and/or  in  a  highly  reverberant  (diffuse)  environment,  the  random  errors 
and  the  phase  mismatch  error  will  become  significant  hence  making  the 
measurement  results  inaccuracy.  However,  providing  the  background  noise  is 
time  stationary,  measurements  can  be  made  to  an  accuracy  of  1  dB  even  when  the 
background  level  exceeds  the  source  level  as  much  as  10  dB  [ref. 6  p.22]. 


2.2.2  Intensity  Mapping 

To  perform  intensity  mapping,  the  surface  area  of  interest  is  first  subdivided  into 
a  grid  and  the  intensity  over  each  point  on  the  grid  is  measured.  The  data 
obtained  is  then  further  processed  so  that  maps  of  intensity  can  be  computed  and 
plotted  out  across  the  entire  grid  for  each  frequency  band  of  interest.  In  this  way 
the  precise  area  of  the  domiiumt  contribution  of  noise  can  be  determined. 

Si^e  intensity  is  a  vector  quantity,  it  is  possible  to  locate  the  sound  sources  and 
sinks  by  mapping  the  source  area;  a  pxjsitive  intensity  indicates  a  source,  and  a 
negative  intensity  indicates  a  sink.  It  is  quite  possible  to  find  a  source  next  to  a 
sink  on  the  same  machine  [ref  .2]. 
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2.2.3  Sound  Transmission  Loss  (Reduction  Index) 

While  Iho  ciinventional  procedure  for  the  measurement  of  sound  transmission 
loss  requires  the  test  specimcnt  to  be  placed  between  two  reverberant  rooms 
(called  a  "transmission  suite").  The  measurement  of  sound  intensity  needs  only 
one  reverberant  room  (rcf.2,7];  in  this  procedure  the  source  is  placed  in  a 
reverberant  room  to  provide  a  diffuse  incident  field,  and  the  transmitted  power  is 
measured  using  sound  intensity  technique. 


2.2.4  Source  Ranking 

A  source  within  a  complicated  structure  may  radiate  sound  in  some  components 
and  absorb  sound  in  others.  In  order  to  effectively  reduce  the  noise  level  of  this 
type  of  source  it  is  necessary  to  rank  the  noise  level  for  each  of  the  source 
components  individually  and  then  focus  on  the  source  component  that  makes  a 
dominant  contribution  to  the  total  noise  level.  This  makes  the  measurement  of 
sound  intensity  a  powerful  tool  because  one  can  define  a  measurement  surface 
that  encloses  the  component  of  interest,  and  treat  all  other  noise-  radiating 
components  as  background  noise  as  long  as  the  noise  is  stationary.  Further,  the 
total  sound  power  can  be  calculated  simply  by  adding  the  partial  sound  powers 
from  all  the  components  that  radiate  noise  [ref .8]. 


3.  Calculation  of  sound  intensity 


In  general,  the  component  of  intensity  in  the  r-direction  is  dehned  as  the  amount 
of  energy  passing  through  an  unit  area  per  unit  of  time; 

I  =  —  (!' 

'At  '  ' 

where  A  and  t  are  the  energy,  the  area  and  the  time  (in  seconds). 

In  acoustics,  energy  is  the  work  done  by  the  sorind  field  on  the  air  particles 
causing  them  to  move: 

E  =  Work  =  Fd, 

where  F  and  dr  are  the  force  and  the  distance, 
since  F/A  is  the  sound  pressure  p,  we  have: 

Er  =  p.A.dr  (2) 

Therefore  to  describe  intensity  (2)  is  substituted  into  (1) : 

pAd,  pd, 

'  At  t 
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Ir  =  PU, 


0) 


since  iir  =  dr  / 1  is  the  air  particle  velocity  in  the  r-direction. 

Equation  (3)  above  is  the  expression  for  an  instantaneous  intensity  at  a  point  r  in 
space.  For  the  mean  intensity  we  have: 

<Jr>  =  P“,  (4) 

where  the  indicates  time  averaging. 

In  three  dimensions,  the  Intensity  vector  can  be  expressed  as: 

<l>  =  pu  (5) 


4.  The  Two  Microphone  Method  For 
Practical  Sound  Intensity  Measurements 


Equation  (5)  shows  that  sound  intensity  is  the  time  averaged  product  of  the  sound 
pressiire  and  air  particle  velocity  at  a  point  in  space.  In  principle,  a  sound 
intensity  measuring  device  should  consist  of  transducers  for  the  detection  of  the 
air  particle  velocity  and  the  sound  pressure  signals  (these  transducers  are 
associated  in  an  intensity  probe).  The  signals  from  the  transducers  are  then 
multiplied  and  time  averaged  to  determine  the  mean  sound  intensity  in  the 
direction  of  the  probe  axis. 

Currently,  two  methods  are  widely  used  for  practical  measurements  of  sound 
intensity.  The  first  method  uses  a  particle  transducer  and  a  sound  pressure 
transducer  (the  p-u  probe).  The  second  method  employs  two  soimd  pressure 
transducers  (the  p-p  probe).  In  the  latter  method,  the  particle  velocity  is 
determined  from  the  spatial  pressure  gradient  via  the  Euler  equation,  using  a 
finite  difference  approximation,  while  the  pressure  is  approximated  as  the 
average  of  the  two  transducer  pressures.  The  approximations  involved  in  the  p-p 
probe  are  discussed  below. 


4.1  The  Estimator  for  Particle  Velocity  (tif) 

According  to  Newton's  second  law  of  motion: 


F  =  ma 


(6) 


Applying  this  law  to  a  unit  volume  of  air  yields  the  Euler  equation: 


P 


du 

dt 


(6a) 


Integration  (6a)  with  respect  to  time  will  give  the  air  particle  velocity: 


u 


(7) 


where  p  is  the  density  of  the  air,  u  is  the  air  particle  velocity  and  Vp  is  the 
pressure  gradient  at  a  point  in  space.  Then  from  (7)  the  component  of  u  in  the 
r  direction  is  given  by: 


u 


r 


-1/p 


J 


(  dp/ dr  )  dt 


(8) 


By  using  the  finite  difference  approximation  method,  the  pressure  gradient  dp/ dr 
can  be  estimated  in  practice  by  measuring  the  pressures  p^  and  p|,  at  two  closely 
spaced  points  separated  by  a  distance  Ar : 


^  -  Pfa  •  P. 

dt  Ar 


(9) 


Note  that  this  approximation  is  valid  only  if  Ar  is  small  compared  with  the 
shortest  wave  lengths  in  the  measured  sound  field. 

Substituting  equation  (9)  into  equation  (8)  the  estimator  fi,.  for  the  particle 
velocity  [ref.5]  is; 

Ur=^J(Pb-P.)dt  (10) 


4.2  The  Estimator  for  Sound  Pressure  (P,) 


Using  a  system  of  two  microphones  for  estimation  of  above,  the  sound 
pressure  p  can  be  estimated  as  the  average  of  the  pressure  pg  and  pj,.  Hence: 


p  Pb+P, 
'■  2 


(11) 
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4.3  The  Estimator  for  Sound  Intensity  (ty) 

Finally,  the  estimated  soimd  intensity  in  the  r-direction ,  !f  is  given  by: 

(Pb-PJdt  (12) 

zpAr  * 

Note  that  equation  (12)  is  used  in  the  analog  signal  processing  analyser  in  the 
calctilation  of  intensity. 

In  the  frequency  domain,  the  estimated  sound  intensity  can  be  calculated 
from  the  imaginary  part  of  the  cross  power  spectrum,  of  the  signals  from  the 
two  microphones  [ref.l] 

This  applies  to  an  ideal  continuous  frequency  spectrum.  In  practical  measurement 
and  aiudysis  a  luurow-band  frequeiKy,  dual  channel  FFT  anal3rser  can  be  used  for 
calculating  the  intensity.  In  such  cases,  the  equation  (13)  becomes: 


I  _  -1  ^Iin[G^(nAf)]  .  . 

'“2itfpAr;^  nAf  '  ' 

where  N  is  the  number  of  spectral  lines  in  the  cross  power  spectrum  and  Af  is 
the  frequency  increment  (resolution)  between  the  spectral  lines. 


5.  Sound  Intensity  Measurement  System 


A  practical  system  for  sound  intensity  measurement,  using  the  frequency  domain 
analysis  derived  above,  is  described  below.  The  system  comprises: 

-  P-P  type  sound  intensity  probe  (  Bruel  &  I^er  type  3520 ). 

-  A  phase  matched  microphoite  pair  (  Bruel  &  Kjaer  type  4183 ). 

-  Dual  chaimel  FFT  analyser  ( Ono  Sokki  CF-350 ). 

-  A  microcomputer  for  post  processing  of  data. 

The  sound  intensity  probe  (t)rpe  3520)  operates  with  a  pair  of  microphones  (type 
4183)  separated  by  a  solid  spacer  (see  Appendix  1(A)).  The  4183  consists  of  tivo 
phase-matched  and  amplitude  matched,  prepolarised  microphones  which  feature 
special  phase-corrector  units.  These  microphones  are  used  to  measure  the  sound 
pressures  Pa  and  Pb  at  two  closely  spaced  points  and  the  spacer  provides  Ar,  the 
separation  distance  between  these  two  points. 

The  Ono  Sokki  CF-350  FFT  analyser  is  used  to  convert  the  time  domain  into  the 
frequency  domain  by  taking  a  Fourier  Transform  of  the  output  signals  from  the 


L,., 


13 


two  microphones.  It  then  produces  and  stores  the  cross  spectrum  of  the  two 
signals  for  each  point  of  measurement. 

A  microcomputer  is  used  for  interfacing  and  processing  the  cross  spectrum  data 
obtained  from  the  Ono  Sokki  CF-350. 

A  schematic  diagram  of  this  sound  intensity  measurement  system  is  presented 
in  Figure  1.  The  photographs  of  the  complete  system  are  shown  in  Appendix  1(C). 
The  limitations  of  this  system  are  discussed  in  section  9. 


Figure  1:  Schematic  diagram  of  a  sound  intensity  measurement  system. 


6.  Computer  Program  For  Calculation  Of 
Sound  Intensity  From  Cf-350  Fft  Analyser  Data 


From  equation  (14),  the  sound  intensity  at  a  point  r  in  space  can  be  deduced  from 
the  cross  spectrum  of  the  two  microphones  (intensity  probe)  using  a  dual  channel 
FFT  analyser. 

A  computer  program  has  been  developed  for  processing  sound  intensity 
measurement  data  (in  the  form  of  cross  spectrum)  using  the  programming 
language  ASYST.  The  program  performs  the  following  functions: 

(1)  It  interfaces  with  the  Ono  Sokki  (CF-350)  dual  channel  FFT  analyser  via  a 
GPIB  board.  This  enables  the  reading  of  cross  spectrum  data  from  the 
analyser  block  memory  into  the  PC.  Note:  each  cross  spectrum 
corresponds  to  a  measurement  point  on  the  control  surface. 

(2)  It  calculates  the  sound  intensity  from  the  cross  spectrum.  The  calculation  is 
based  on  equation  (14)  of  section  4.3  and  the  computed  intensity  represents 
the  component  which  is  normal  to  the  measurement  surface  (also  the 
direction  of  the  probe  axis). 

(3)  It  provides  post  processing  of  sound  intensity  data  such  as: 

-  Calculation  of  total  radiated  power  by  integrating  of  partial  sound 
powers  (from  sub-areas  of  the  control  surface). 

-  3-dimensional  plotting  of  intensity  distribution  over  the 
area  of  the  measurement  surface,  and 

-  Third  octave  band  analysis  of  the  power  radiated  from  the 
control  surface. 

Note:  the  formation  of  third  octave  bands  from  narrow  bands 
(spectrum  lines)  includes  Hanning  window  compensation. 

The  method  used  is  based  on  that  of  the  CF-350  analyser. 

(4)  It  allows  the  user  to  re-execute  the  program  for  another  set  of  data. 

(5)  It  provides  print-out  functions  for  all  plots  and  result  listings. 

Some  of  the  program  features  are: 

•  menu  selection  for  driving  the  program  as  shown  in  Figure  2, 

•  screens  for  data  entry.  Usually  there  is  a  brief  description  about  the  process 
and  types  of  data  to  be  input  (Fig.3), 

•  status  of  current  process,  instructions  for  loading  of  data  etc.  (Fig.4), 

•  in  third  octave  band  analysis,  the  program  displays  the  centre  frequency, 
band  number,  and  band  value  both  in  engineering  units  and  in  dB  (with 
reference  to  lOE-1 2  W  )  (Fig.5), 
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a  3-dunensional  plot  of  sound  intensity  distribution  (intensity  map)  over  a 
quasi-surface  area  as  shown  in  Figure  6, 

a  sample  histogram  plot  of  analysis  results  in  third  octave  band,  as  is  shown 
in  Figure  7. 
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Figure  2:  Menu  selection  for  driving  the  sound  intensity  jrrogram. 
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Figure  6:  A  3-dimensional  plot  of  sound  intensity  distribution  (intensity  map)  oi>er  the 

measured  surface. 
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Figure  7:  A  histogram  plot  from  a  third  octave  band  analysis  results. 
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7.  Errors  In  The  Estimation  Of  Sound 
Intensity  With  P-P  Type  Intensity  Probe 


7.1  Finite  Difference  Approximation  Error 

As  described  in  section  4,  to  estimate  the  air  particle  velocity,  the  pressure 
gradient  dp/di  has  been  approximated  by  Ap/Ar  using  a  finite  difference 
method.  This  approximation  is  valid  when  Ar  is  small  compared  to  the  shortest 
wavelength  in  the  sound  field  to  be  measured.  At  high  frequency,  when  the 
corresponding  wavelength  is  small  compared  to  the  effective  microphone 
separation  then  dp /dr  *  Ap/Ar  and  the  estimate  of  intensity,  I,  will  become 
inaccu  irate. 

For  a  particular  microphone  separation,  Ar,  there  will  be  an  upper  hrequency 
limit  of  the  measurement  system  beyond  which  results  may  be  inaccurate.  As  an 
example,  the  upper  frequency  limits  for  the  standard  microphone  separations  of 
the  Bruel  &  Kjaer  sound  intensity  probe  beyond  which  the  error  is  greater  than 
1  dB  are  ( assuming  plane  wave  along  the  probe  axis  [ref.  9  pp.l3]) : 


Ar  fmmi 


fmax.  (VHx) 


6 

12 

50 


10 

5 

1.2 


From  the  table  above  we  can  see  that  if  a  smaller  spacer  is  chosen,  the  upper 
frequency  limit  of  the  system  can  be  increased.  For  an  error  within  1  dB,  the 
smallest  wavelength  measured  in  the  sound  field  should  be  at  least  6  times  the 
microphone  spiacer  Ar. 


7.2  Probe  Diffraction  Error 

Due  to  the  presence  of  the  probe,  the  sound  held  will  be  distorted  by  diffraction 
and  there  will  be  variations  in  the  sensitivity  and  acoustical  separation  of  the 
microphones  [ref. 9  ppl2].  This  places  an  upper  frequency  limit  on  practical 
probes.  Usually  the  evaluation  of  probe  diffraction  characteristics  is  carried  out  by 
the  manufacturer  and  the  probe  frequency  limit  is  stated  with  the  hardware. 
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7.3  System  Phase  Mismatch  Error 
7.3.1  System  Phase  Mismatch 


It  can  be  shown  that  the  intensity  measured  at  a  point  in  the  sound  field  is  directly 
related  to  the  phase  difference  <I>gb  detected  ( by  the  system )  between  the  two 
microphones  of  the  sound  intensity  probe.  Ideally,  this  phase  angle  should  purely 
be  the  phase  change  of  the  sound  field  pressure  across  the  two  microphones  of  the 
probe.  In  practice,  there  always  exists  a  phase  mismatch  in  all  sound  intensity 
measurement  systems.  This  system  phase  mismatch  is  the  combined  effect  of: 

•  phase  mismatch  between  the  two  microphones  of  the  probe, 

•  phase  mismatch  between  chaimels  of  the  analyser. 

Hence  in  an  actual  measurement,  the  phase  difference  detected  between  the 
two  microphones  is  the  sum  of  the  actual  phase  change  of  the  sound  field  and  the 
system  phase  mismatch. 

For  a  measurement  to  be  accurate,  the  system  phase  mismatch  must  be  kept 
small.  This  can  be  achieved  by  using  an  analyser  with  two  highly  phase  matched 
chaimels.  Also  the  probe  should  only  employ  a  phase  matched  microphone  pair, 
such  as  those  especially  designed  for  sound  intensity  measurements. 


7.3.2  Phase  Mismatch  and  Error  at  Low  Frequency 

For  the  error  due  to  phase  mismatch,  in  the  estimate  of  intensity  1,  to  be  negligible 
the  phase  change  of  sound  pressure  across  the  microphones  must  be  many  times 
larger  than  the  system  phase  mismatch.  This  is  analogous  to  the  signal  to  noise 
ratio  in  an  electrical  system.  Consequently  the  effect  of  system  phase  mismatch  is 
most  critical  for  small  microphone  spacings  and  at  low  frequency  since  the  sound 
field  phase  change  is  small  in  these  cases  (ref.4  p.ll5]. 

To  reduce  the  effect  of  phase  mismatch  at  low  firequency  a  larger  spacer  can  be 
used,  but  this  reduces  the  system  upper  frequency  limit,  as  shown  in  section  7.1. 
Therefore  at  low  frequency  a  large  spacer  should  be  used  and  at  higlh  frequency  a 
small  spacer  is  preferred. 

Hence  in  a  sound  intensity  measurement  system,  for  a  chosen  Ar  there  is  a  low 
frequeiKy  limit  beyond  which  the  error  due  to  system  phase  mismatch  is 
uiuicceptable  in  the  estimate  of  I. 


7,3.3  Phase  Mismatch  and  The  Pressure-Residual  Intensity  Index 

If  a  sound  wave  is  incident  at  90°  to  the  probe  axis,  the  two  microphones  are 
exposed  to  the  same  sound  signal .  In  this  case  the  field  pressure  phase  difference 
across  the  microphones  is  zero  and  any  phase  difference  detected  is  the  system 
phase  mismatch.  The  intensity  corresponding  to  this  phase  mismatch  is  called  the 
residual  intensity.  The  residual  intensity  depends  on  both  the  magnitude  of  the 
system  phase  mismatch  and  the  sound  pressure  at  the  microphones.  It  can  be 
shown  that  for  a  chosen  microphone  spacing  Ar  and  at  a  given  frequency,  f ,  the 


difference  between  pressure  level  and  residual  intensity  level  is  unchanged.  This 
level  difference  is  defined  as  the  pressure-residual  intensity  index  Spig- 

SpIO  =  Lp(90)  -  Lires(90)  05) 

where  Lp^9Q^  ^lres(90)  Ih®  sound  pressure  level  and  its  corresponding 

residual  intensity  level  when  the  sound  wave  is  90°  to  the  probe  axis. 

The  pressure-residual  intensity  index  describes  the  phase  mismatch 
characteristics  of  a  particular  measurement  system. 

Note:  One  way  to  measure  the  pressure-residual  intensity  index  SpjQ  is  to  use  an 
acoustic  coupler  where  a  plane  wave  incident  at  90°  to  the  probe  axis  can  be 
simulated. 


7.4  Air  Flow  Disturbance 

The  disturbance  caused  by  an  air  flow  can  contaminate  the  signals  horn  sound 
intensity  probes  [ref.  4  p.l24].  Therefore  windscreens  should  always  be  used  for 
an  outdoor  measurement  or  where  there  is  an  .dr  flow  within  the  vicinity  of  the 
measurement  surface. 


7.5  Random  Error 

If  the  sound  field  is  contaminated  with  extraneous  noise  sourcefs)  aivu/or  high 
diffuse  background  noise  the  random  error  in  the  estimate  intensity  can  be  severe. 
It  has  been  shown  by  Gade  [ref.  12]  that  in  a  partially  diffuse  sound  field  the 
random  error  depends  on  the  BT-product  and  the  field  measured  pressure- 
intensity  index,  Opj.  For  a  measurement,  if  the  value  of  5pj  is  large  then  the 
averaging  time  T  must  increase  in  order  to  reduce  the  random  error  in  the 
measured  intensity. 

Note:  BT-product  is  the  product  of  the  frequency  bandwith  B  and  the  averaging 
time  T. 


7.6  Other  Errors 

If  the  output  of  the  sound  source  is  not  stationary  with  time  or  extraneous  noises 
are  transient  then  there  will  be  an  error  in  the  intensity  measured. 
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8.  Pressure-Intensity  Index  And  The  System 
Dynamic  Capability 


8.1  Pressure-Intensity  Index 

The  pressure-intensity  index  is  defined  as: 

Spl  =  ^)  -  M  (^6) 

where  Lp  is  the  pressure  level,  and  L|  is  the  measured  intensity  level  at  the  point 
of  measurement. 

Sp{  is  a  measure  of  the  ratio  between  the  true  hee  held  intensity  I  to  the  measured 
intensity  I,  in  dB  [ref.  4  p.117].  Therefore  6p|  should  be  as  small  as  possible  so 
that ! » I.  A  small  measured  intensity !  will  correspond  with  a  large  value  of  5p[. 

If  !  is  small  enough,  the  error  due  to  system  phase  mismatch  will  become 
significant,  hence  making  the  measurement  inaccurate.The  effect  of  phase 
mismatch  error  on  the  measured  intensity  is  determined  by  the  pressure  intensity 
index  Spj  and  the  system  dynamic  capability  Lj  (to  be  discussed  later).  Beside 
this,  the  random  error  in  the  estimate  intensity  is  also  dependent  on  the  [ref.  4 
pp.138,140 ;  ref.l2  p.l5].  A  large  value  of  8pj  will  correspond  with  a  high  level  of 
difficulty  in  making  an  accurate  measurement  of  sound  intensity. 

The  6p]  can  be  reduced  by  placing  the  probe  closer  to  the  source  to  improve  the 
signal  to  noise  ratio  or  placing  sound  absorption  materials  around  the  walls 
(outside  the  measurement  surface)  to  reduce  the  reflections  of  sound  waves  at  the 
boundaries  etc. 


8.2  System  Dynamic  Capability 

For  the  measured  intensity  to  have  a  reasonable  level  of  accuracy,  the  actual  phase 
change  of  the  sotmd  field  across  the  microphones  must  be  large  compared  to  the 
system  phase  mismatch.  This  is  equivalent  to  the  pressure-residual  intensity  index 
5p|o  being  much  larger  than  the  pressure-intensity  index  Spj.  For  an  accuracy  in 
the  measurement  of  intensity  to  within  1  dB  and  0.5  dB  the  held  measured  5pj 
must  be  smaller  than  (8pjg  -  7)  dB  and  (6p{o  - 10)  dB  respectively.  From  this,  the 
system  dynamic  capability  is  defined  as: 

Ld  “  Spio  -  K  (17) 

where  K  is  a  constant  which  is  dependent  on  the  level  of  accuracy  to  be  achieved 
(  eg.  7  dB,  10  dB  for  an  accuracy  of  ±  1  dB,  ±  0.5  dB  respectively) . 

The  field  measured  Spj  miist  not  exceed  the  level  indicated  by  in  order  to 
achieve  the  level  of  accuracy  proposed  by  the  constant  K.  Usually  the  phase 
mismatch  error  is  significant  at  low  frequency.  The  fi:equency  in  which 
Lj  ^  5pi  is  regarded  as  the  low  frequency  limit  of  the  system. 
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9.  Eva  lua  tion  Of  System  Performance 


Following  is  a  description  of  the  limits  and  capabilities  of  the  system 
configuration  as  described  in  section  5  (Sound  intensity  measurement  system). 

(1)  High  frequency  limit 

The  Bruel  &  Kjaer  intensity  probe  type  3520  together  with  the  phase  matched 
microphone  pair  type  4183  ( 12  mm  and  50  mm  spacers )  has  an  upper  frequeiKy 
limit  of  5  kHz  (ref.  14].  This  determines  the  5  kHz  upper  h^uency  limit  for  the 
measurement  system  as  a  whole. 


(2)  The  Processor  Real  Time  Analysis 

The  Ono  Sokki  dual  chaimel  FFT  analyser  (CF-350)  can  operate  at  a  frequency 
range  up  to  40  kHz.  For  real  time  analysis,  2  kHz  is  the  range  limit. 


(3)  The  Processor  (CF-350)  Channel  Phase  Mismatch 

By  feeding  common  electrical  signals  to  the  CF-350  channel  inputs,  the  phase 
mismatch  between  the  two  channels  of  the  CF-350  has  been  evaluated  using  the 
phase  part  of  the  cross  spectrum  ( or  transfer  function).  It  was  found  that  for  a 
5  kHz  frequency  range  and  a  frequency  resolution  of  12.5  Hz  (400  line  spectrum) 
the  CF-350  has  a  typical  phase  mismatch  which  was  less  than  or  equal  to  ±  0.2° 
between  its  two  channels  (random  signal  inputs  with  1024  times  of  averaging). 


(4)  The  Microphone  pair  (4183)  Phase  Mismatch 

Typically  the  type  4183  has  a  phase  match  which  is  better  than  0.2°  from  40  Hz  to 
700  Hz  and  it  is  better  than  (f  /  3500)°  for  frequency  f  up  to  5  kHz  [ref.  14].  The 
calibration  of  the  phase  part  (supplied  by  the  manufacturer)  for  the  microphone 
pair  which  is  used  in  our  system  is  shown  in  the  Appendix  1(B). 


(5)  System  Pressure-Residual  Intensity  Index  and  Dynamic  Capability 

Once  soimd  intensity  calibration  has  been  completed,  the  type  3541  Bruel  &  Kjaer 
sound  intensity  calibrator  can  be  used  for  measurement  of  the  system  pressure- 
residual  intensity  index.  In  the  acoustic  coupler  chamber  (UA  0914),  the  two 
microphones  of  the  probe  are  exposed  to  the  same  source  of  pink  noise  (ZI  0(B5). 
The  pressure-residual  intensity  index  can  be  computed  by  subtracting  the 
detected  residual  intensity  level  from  the  sound  pressure  level  at  the  two 
microphones. 

Figure  8  gives  typical  pressure-residual  intensity  indices  of  the  system  with  the 
microphone  pair  t)rpe  4183,  a  12  mm  microphone  spacer,  and  a  frequency  range  of 


0-5  kHz.  The  dynamic  capability  of  the  system  with  K  =  7  dB  and  10  dB  are  also 
presented  in  this  figure. 
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System  Res.  P-l  index  &  Dynamic  capabiiity 
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Figure  8:  The  SpiQ  ^  and  L^for  K=  10,  7  dB  offhc  system, 
(a)  Numerical  values,  (b)  Graphical  presentation. 


^Note:  since  the  option  for  calculation  of  pressure-residual  intensity  index  is  not 
yet  available  with  the  oirrent  version  of  sound  intensity  program  the  mean 
pressure  level  has  been  obtained  manually  hvm  the  CF-350  then  substracted  by 
the  residual  intensity  level  which  is  calculated  firom  the  soimd  intensity  program. 
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10.  Procedures  For  Calibration  Of  The  System 


For  reliable  results,  a  sound  intensity  measurement  system  should  be  calibrated 
properly  before  any  measurements  are  made.  The  intensity  can  be  calibrated  in  an 
anechoic  chamber,  a  plane  wave  tube  or  an  acoustic  coupler.  At  MRL  the 
measurement  system  is  calibrated  with  the  Bruel  &  Kjaer  sotind  intensity 
calibrator  type  3541.  This  type  of  calibration  is  carried  out  in  an  acoustic  coupler 
where  sound  waves  of  both  0  and  90  incidence  to  the  probe  axis  can  be 
simulated.  These  allow  us  to  calibrate  both  pressure  and  sound  intensity 
sensitivities  of  the  system  and  also  allow  the  measurement  of  system  pressure- 
residual  intensity  index,  5p]Q. 

Following  is  an  outline  of  the  procedure  for  calibration  of  the  system,  with 
respect  to  the  calibrator  type  3541  mentioned  above.  For  further  details  such  as 
the  calculations  of  the  correction  terms  for  the  actual  ambient  conditions, 
operation  of  the  CF-350  etc.,  the  reader  is  referred  to  the  manuals  of  the 
corresponding  instruments. 


10.1  Pressure  Calibration 

1.  Calculate  the  correct  pressure  value  given  by  the  calibrator  under  the  actual 
ambient  conditions  of  the  measurement.  This  is  done  by  applying  the 
correction  terms  (specified  in  the  instruction  manual)  to  the  reference 
pressure  stated  in  the  calibration  chart  of  the  pistonphone. 

2.  Set  the  CF-350  to  operate  at  500  Hz  frequency  range  and  2  volts  amplitude 
range  for  both  channels. 

3.  Set  the  CF-350  to  operate  on  third  octave  band  mode. 

4.  Insert  the  microphones  into  the  coupler  ports  which  are  intended  for  pressure 
calibration.  This  enables  the  microphone  pair  to  be  exposed  to  the  same  soimd 
signal  from  the  source. 

5.  Place  the  pistonphone  type  4228  on  the  coupler  and  turn  it  on.  This  gives  a 
reference  signal  of  X  Pa  (Y  dB  ref  ^  p  Pa  )  at  251.25  Hz  inside  the  coupler 
chamber.  Where  X  and  Y  are  the  corrected  values  of  the  reference  signal 
calculated  in  step  1. 

6.  Set  the  CF-350  to  display  the  power  spectrum  of  channel  A  and  move  the 
cursor  to  the  250  Hz  peak. 

7.  Use  the  soft  key  [SP/EU]  of  the  CF-350  to  set  the  sensitivity  for  channel  A  so 
that  the  current  cursor  corresponds  to  X  Pa.  A  value  will  be  assigned  to  the 
EU/V  on  the  screen,  this  is  the  sensitivity  factor  for  channel  A. 

8.  Set  the  CF-350  to  display  the  power  spectrum  of  channel  B,  move  the  cursor 
to  the  250  Hz  peak  then  set  the  channel  B  sensitivity  similarly  to  step  7  above. 

10.2  Intensity  Calibration 

1.  For  the  particular  microphone  spacing  at  the  probe,  calculate  the  correct 
intensity  level  given  by  the  calibrator,  under  the  actual  ambient  conditions. 
Here,  correction  terms  (specified  in  the  instruction  manual)  are  applied  to 
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give  the  values  appropriate  to  the  ambient  conditions  rather  than  reference 
conditions  assumed  (stated  in  the  cahbration  chart  of  the  pistonphone). 

2.  Set  the  CF'350  to  operate  at  500  Hz  frequency  range  and  2  volts  amplitude 
range  for  both  channels. 

3.  Set  the  CF-350  to  display  the  cross  power  spectrum  . 

4.  Insert  the  microphones  into  the  coupler  ports  which  are  intended  for  intensity 
calibration.  This  enables  the  sound  source  to  simulate  a  plane  wave  which  is 
0  incident  on  the  probe  axis. 

5.  Place  the  pistonphone  type  4228  on  the  coupler  and  turn  it  on. 

6.  Run  the  program  for  soiuid  intensity  calculation. 

7.  Choose  third  octave  analysis  option  from  the  program  main  menu  to  display 
the  ccdculated  data  in  third  octave  band. 

8.  Check  if  the  intensity  measured  by  the  system  at  the  250  Hz  band  is  matched 
with  the  corrected  intensity  calculated  in  step  1. 


10.3  Measurement  of  Pressure-Residual  Intensity  Index  SpjQ 

1.  Set  the  CF-350  to  operate  on  frequency  range  of  interest  and  adjust 
amplitude  range  of  both  channels  to  give  the  optimum  signal  to  noise  ratio. 

2.  Set  the  CF-350  to  display  the  cross  power  spectnun  and  set  the  number  of 
averaging  N  («  512  - 1024). 

3.  Insert  the  microphones  into  the  coupler  ports  which  are  intended  for  sound 
pressure  calibration.  This  enables  the  sound  source  to  simulate  a  plane  wave 
which  is  90  incident  on  the  probe  axis. 

4.  Place  the  source  of  pink  noise  ZI 0055  on  the  co.  pler  and  turn  it  on. 

5.  Start  the  averaging  process  for  measurement  ;  toss  spectrum  across  the 
microphones  until  it  finishes. 

6.  Rxin  the  program  for  sound  intensity  calculation.  When  asked,  input  the 
microphone  spacing  to  be  used  with  the  probe  and  aU  other  reqtiired 
parameters. 

7.  Choose  third  octave  analysis  option  from  the  program  main  menu  to 
display  the  calculated  data  in  third  octave  band.  This  is  the  residual 
intensity  spectrum  of  the  system  with  respect  to  the  pressure  produced  by 
the  broad  band  sound  source  ZI  0055. 

8.  Display  either  power  spectrum  of  channel  A  or  B  in  third  octave  band.  The 
reading  dB  value  on  the  CF-350  is  referenced  to  IPa.  Add  94  dB  to  this  value 
to  obtain  the  sound  pressure  level  with  reference  to  20  ^  Pa.  For  each 
frequency  band,  subtract  the  residual  intensity  level  from  the  pressure  level. 
This  level  difference  between  sound  pressure  level  and  residual  intensity  is 
the  pressure-residual  intensity  index  Spio  ■ 


11.  Sound  Power  Determination  Of  A 
Reference  Source 

To  show  the  application  of  sound  intensity  measurement  in  sound  power 
determination,  the  sound  power  output  from  a  reference  source  (Bruel  &  Kjaer 
type  4205)  has  been  measured  using  the  intensity  measurement  system  in  the 
configuration  which  was  described  in  section  5  (Sound  intensity  measurement 
system). 

The  power  output  of  the  sound  source  type  4205  can  be  varied  continuously 
between  40  and  approximately  100  dB  re  1  pW.  The  output  level  can  be  broad 
band  pink  noise  from  100  Hz  to  10  kHz  range  or  octave  band  filtered  noise  by 
using  one  of  the  7  built-in  octave  band  pass  filters.  Because  the  upper  frequency 
limit  of  the  system  is  5  kHz,  we  cannot  use  the  reference  output  given  by  the 
broad  band  pink  noise.  Instead,  the  octave  band  pass  filters  (from  125  Hz  to  4  kHz 
centre  frequency)  were  used  to  give  a  nominated  output  of  85  dB  re  1  pW  from 
each  band. 

(1)  Equipment  used 

•  The  intensity  measurement  system  consists  of  a  microcomputer  to  run  the 
sound  intensity  program,  the  CF-350  dual  channel  FFT  analyser,  a  Bruel  & 
Kjaer  sound  intensity  probe  type  3520,  a  phase  matched  microphone  pair  type 
4183  with  a  12  mm  or  50  mm  spacer. 

•  The  sound  power  source  type  4205  consists  of  two  separate  units:  the 
generator,  containing  all  the  controls,  filters,  amplifiers,  level  meter  etc.  and 
the  sound  source  HP  1001  containing  two  loud  speakers  with  the  associated 
crossover  networks. 


(2)  Calibration  of  the  equipments 

•  The  source  type  4205  has  been  calibrated  according  to  the  procedures 
described  in  its  instruction  manual.  Due  to  the  equipment  availability  the 
source  type  4205  has  been  calibrated  using  the  Bruel  &  Kjaer  real  time 
frequency  analyser  type  2144  and  one  of  the  microphone  from  the  probe. 

•  The  sotmd  intensity  measurement  system  has  been  calibrated  for  both  sotmd 
pressure  and  intensity  sensitivities  according  to  the  procedures  described  in 
section  10  (system  calibration  procedure)  using  the  Bruel  &  Kjaer  sound 
intensity  calibrator  type  3541. 

i 

(3)  The  measurement  surface  and  tiie  environments 

The  measurement  has  been  carried  out  in  a  room  of  8.0  x  8.0  x  3.6  m 
approximately.  The  floor  and  the  ceiling  are  rigid  concrete,  and  the  floor  is  tiled 
with  vinyl  sheet.  The  room  also  contains  some  timber  cabinets,  equipments  etc. 

There  was  a  low  level  of  background  noise  during  the  measurements.  This  was 
the  fan  noise  fiom  the  corridor  outside. 
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The  measurement  (control)  surface  is  an  imaginary  cube  of  dimensions 
1  X  1  X  1  m  .  The  source  was  placed  on  the  floor  at  the  center  of  the  cube's  bottom 
face.  The  sotmd  intensity  components  normal  to  its  faces  was  to  be  measured. 
Data  were  taken  on  all  of  the  cube's  faces  except  the  bottom  one  since  it  is 
assumed  to  be  reflecting  sound  energy  back  to  the  volume  enclosed  by  the  cube 
surface. 

Each  face  of  the  cube  was  subdivided  into  a  grid  of  4  (2  x  2)  elements  and  each 
had  an  area  of  0.25  m^  (0.5  x  0.5  m).  This  makes  a  total  of  20  measiuement  points 
over  the  whole  cube.  The  normal  component  of  sound  intensity  at  the  center  of 
each  element  was  measured  and  the  power  radiated  from  each  of  these  elements 
is  given  by  the  product  of  its  area  and  the  normal  intensity  component.  The  total 
sound  power  radiated  from  the  cube  can  be  calculated  by  integration  of  these 
elemental  sound  powers. 


(4)  Measurement  settings 

Settings  of  the  sound  source: 

•  For  each  octave  band  from  125  Hz  to  4  kHz ,  the  sound  source  was  set  to  give 
an  output  power  level  of  85  dB  with  reference  to  1  pW. 

Settings  of  the  FFT  analyser: 

•  The  engineering  unit  calibration  factors  for  the  two  channels  were  set  during 
the  calibration  process. 

•  For  each  measurement  band,  the  frequency  range  on  the  CF-350  was  set  so 
that  it  gave  the  highest  frequency  resolution  on  the  cross  spectrum  between 
the  two  channels.  In  this  way  we  obtained  the  most  number  of  spectrum  lines 
within  the  band  of  interest  and  consequently  the  result  was  more  accurate 
when  these  lines  were  grouped  to  synthesise  the  corresponding  octave  band 
value. 

•  The  Hanning  window  was  used  to  reduce  the  leakage  effects  in  the 
measurement  data.  During  synthesis  of  the  third  octave  bands  the  effect  of 
this  applied  window  was  compensated  for  by  applying  the  0.66  factor 
described  in  the  CF-350  manual. 

•  The  amplitude  range  of  the  cross  spectrum  was  set  manually  until  an 
optimum  signal  to  noise  ratio  obtained.  The  number  of  averages  was  set  to 
256. 


(5)  Discussion  of  results 

The  total  sound  power  radiated  from  the  cube  was  computed  from  the  intensity 
measurements  over  the  cube  (control)  surface  as  described  in  (3)  above. 

The  residts  of  the  measurement  (with  a  12  mm  microphone  spacer)  are 
summarised  in  Table  1  and  Figure  9.  More  details  about  the  measurements  can  be 
found  in  Appendix  2. 
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Figure  9:  Sound  power  measurement  of  the  refn’ence  source  with  85  dB  nominated  for 
each  octave  band  (loith  a  12  mm  spacer). 
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From  Table  1  we  can  see  that  in  general,  except  at  125  Hz  and  4  kHz  bands,  the 
measured  power  level  is  quite  close  and  is  less  than  1 0.5  dB  From  the  nominated 
power  level  (85  dB  re  1  pW)  .  The  difference  between  the  measured  level  and  the 
nominated  level  is  due  to  the  combined  effects  of : 

•  Uncertainty  in  the  output  level  of  the  reference  source.  For  the  output  range 
from  50-90  dB,  according  to  the  source  instruction  manual,  the  accuracy  in 
the  power  level  output  of  the  source  is  ±  1  dB  for  octave  bands  from  250  Hz 
to  2  kHz  and  ±  15  dB  for  octave  bands  outside  this  frequeiury  range.  The 
uncertainty  in  the  output  level  of  the  source  was  shown  in  Figure  9  as  the 
output  upper  and  lower  bound  curves. 

•  Errors  in  the  estimation  of  sound  intensity  using  the  sound  intensity 
measurement  system.  These  errors  are  due  to  the  system  phase  mismatch, 
random  noise,  insuffidence  number  of  measurement  points  in  the  estimation 
of  sound  power  output  etc. 

At  4  kHz  octave  band,  the  measured  level  is  0.9  dB  below  the  nominated  value. 
This  is  due  to: 

•  A  large  uncertainty  in  the  source  output,  as  mentioned  above,  for  this 
frequency  band  ( 11.5  dB). 

•  The  upper  frequency  limit  for  4  kHz  octave  band  is  5650  Hz  whereas  the  high 
frequency  limit  of  the  measurement  system  is  5000  Hz.  This  leads  to  the 
underestiination  of  the  mean  preraure  and/or  the  pressure  gradient  and 
consequently  the  intensity  is  also  underestimated  [ref.  1,  9]. 

The  main  reasons  for  the  difference  between  the  measured  level  and  the 
nominated  level  being  greatest  (Z1  dB)  at  the  125  Hz  octave  band  are: 

•  The  uncertainty  in  the  source  output  is  large  for  this  frequency  band  (±1.5 
dB). 

•  The  field  measured  pressure-intensity  index  for  this  band  is  2.4  dB,  whereas 
the  dynamic  capability  for  this  frequency  band  is  only  0.9  dB^,  thus  in  this 
octave  band  the  field  pressrue-intensity  index  exceeds  the  system  dynamic 
capability^.  As  a  consequence,  the  error  caused  by  the  phase  mismatch  in  the 


^Note:  When  a  comparison  is  made  between  6pi  and  L^j,  it  is  imperative  that  5p] 
and  L(]  originate  from  measurements  made  at  the  same  frequency  range  and  then 
identical  frequency  resolution.  This  is  achieved  by  setting  an  identical  range  at  the 
CF-550  analyser  during  the  determination  of  to  that  used  while  making  the 
octave  band  measurements  (determination  of  5p] ). 


^Note:  If  the  error  level  of  1 1  dB  ( instead  of  1 0.5  dB )  is  acceptable  in  our 
intensity  measurements  then  the  system  dynamic  capability  Lj  corresponding  to 
K  =  7  dB  will  have  a  value  of  3.9  dB  for  the  125  Hz  octave  band.  In  this  case,  the 
pressure-intensity  index  5pj  =  2.4  dB  becomes  well  below  the  Lj  value. 
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measured  intensity  level  can  be  exceeded  ±  0.5  dB  ( these  error  bounds  are 
impbed  by  the  system  dynamic  capability  Lj  with  K»10  dB).  This  could  be 
why  there  is  a  fair  difference  between  the  two  leveb  if  it  did  not  mainly  come 
from  the  uncertainty  of  the  source  output. 

To  reduce  the  effect  of  phase  mismatch  at  low  frequency,  we  repeated  the  125  Hz 
octave  band  measurement  with  a  50  mm  spacer.  In  this  case  the  measured  output 
level,  is  86.3  dB,  the  d3mamic  capability,  Lj,  is  7.1  dB  and  the  pressure- 
intensity  index,  Spj,  is  3.3  dB  (Table  2).  Compared  with  Table  1,  we  can  see  that 
although  the  pressure-intensity  index  increased  from  2.4  dB  to  3.3  dB  it  is  now 
well  below  the  system  dynamic  capability  which  is  7.1  dB.  The  measured  output 
level  for  the  125  octave  band  is  now  86.3  dB  with  the  confidence  that  the  phase 
mismatch  error  is  bounded  by  ±  05  dB.  The  measurement  results  corresponding 
with  both  12nun  and  50  mm  spacers  is  presented  in  Figure  10. 


Microphones  spacer 

Ar  «  SOrnm 

Band  fieq. 

125  Hz 

Lw(dB) 

86.3 

I'd  (lOdB) 

7.1 

Spi 

3.3 

Table  2:  Measurements  cf  the  125  octave  band  with  a  50  mm  ndcrophme  spacer. 

Using  a  larger  microphone  spacer  to  improve  the  measurement  results  at  low 
frequency  is  a  possi^ty.  There  are  other  methods  that  can  also  be  used,  such  as: 

•  Implementing  a  computer  procedure  that  will  compensate  for  the  system 
phase  mismatch.  This  helps  to  increase  the  system  pressure-residual  intensity 
index,  SpIO,  hence  the  dynamic  capability,  Lj.  The  mathematical 
principle  of  this  phase  correction  is  described  in  [ref.  4]. 

•  Choosing  a  smaller  measurement  cube  so  that  the  probe  can  be  placed  closer 
to  the  sound  source.  This  improves  the  signal  to  noise  ratio  of  the 
measurement  and  helps  to  reduce  the  field  measured  pressure-intensity 
index,  8pi . 

•  Reducing  the  reverberant  component  of  the  sound  field  under  measurement 
by  placing  sound  absorpbing  materials  on  the  walls  etc.  This  will  also  help  to 
reduce  the  5p|. 

By  manipulation  of  Lj  and  5p|  via  the  methods  described  above,  it  is  possible  to 
use  the  12  mm  spacer  to  measure  intensity  at  relatively  low  hequency  ( eg.  125  Hz 
or  lower  ),  as  long  as  we  have  5p|  i  L^.  The  error  level  associated  with  L^j  (as 
determined  by  K  in  equation  (17))  should  also  be  kept  in  mind. 
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The  above  demonstrates  the  use  of  the  system  dynamic  capability,  Lj ,  the  field 
measured  pressure-intensity  index,  8p|  in  the  evaluation  of  the  sound  field 
quality,  and  the  acceptable  level  of  errors  in  the  measurement  of  sound  intensity. 
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Figure  10:  Final  measurement  results  using  both  12  mm  and  50  mm  spacers. 
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12.  Discussion 


One  of  the  advantages  of  the  sound  intensity  technique  is  that  the  sound  or  noise 
source  of  interest  can  be  measured  on-site  rather  than  in  an  anechoic  duunber  or 
reverberant  room.  Apart  from  the  fact  that  it  is  more  realistic  to  measure  a  noise 
source  from  its  normal  radiating  environment  the  ability  of  on-site  measurement 
also  helps  to  save  the  time  and  effort  for  the  removal  and  reinstallation  of  the 
source.  Further,  using  the  sound  intensity  technique,  sound  potver  determination 
of  a  noise  source  can  be  done  even  in  the  presence  of  other  radiating  sources  as 
long  as  these  background  noises  are  stationary.  These  have  been  demonstrated  in 
our  measurements  of  the  reference  sound  source  Bruel  &  IQaer  type  4206  in 
section  11  where  the  measurements  were  carried  out  in  a  normal  office  space  and 
under  the  presence  of  a  background  fan  noise. 

Sound  intensity  is  a  vector  quantity,  it  has  both  magnitude  and  direction.  At  the 
point  of  measurement,  a  positive  intensity  indicates  the  outward  flow  of  sound 
energy  from  the  surface  of  consideration  while  a  negative  intensity  indicates  an 
energy  flow  in  the  opposite  direction.  This  enables  the  locations  of  the  source  and 
sink  of  a  radiating  surface  to  be  identified,  as  a  source  is  indicated  by  a  positive 
and  a  sink  by  a  negative  intensity.  The  measurements  of  sound  intensity  over  a 
radiating  surface  also  allows  the  construction  of  the  distribution  of  sound 
intensity  over  this  surface.  With  the  aid  of  these  intensity  maps  it  is  possible  to 
locate  the  area(s)  of  strong  noise  radiation.  Sound  intensity  technique  also  has  the 
ability  to  rank  the  sound  source  according  to  radiated  power.  This  can  be  done  by 
dividing  the  source  into  components,  the  sound  power  of  each  component  is 
determined  individually  and  then  compared  and  ranked  in  order  of  sound  power. 

The  sound  intensity  technique  has  some  disadvantage  and  limits  associated 
with  it.  For  a  measurement  ^tem  which  employs  a  p-p  type  probe,  there  is  an 
inherent  systematic  error  since  the  air  particle  velocity  and  the  sound  pressure 
have  been  approximated  by  the  finite  difference  method.  This  was  discussed  in 
section  7.1  and  this  type  of  error  imposes  an  upper  frequency  limit  to  the  system 
for  a  given  microphone  spacer  of  the  probe.  Measurements  attempted  beyond  this 
limit  will  underestimate  the  true  intensity  [ref.l].  The  error  caused  by  the  phase 
misnuitch  of  the  system  is  important  in  the  sound  intensity  measurement 
technique.  This  type  of  error  is  worst  at  low  frequency  when  the  magnitude  of  the 
phase  mismatch  is  ab^ut  the  order  of  the  actual  sound  field  phase  difference 
across  the  two  microphones  of  the  probe.  Consequently,  there  is  a  frequency  limit 
below  which  the  error  in  the  measurement  of  sound  intensity  is  not  acceptable. 
This  low  frequency  limit  is  determined  by  the  pressure-intensity  index  and  the 
system  dynamic  capability  as  discussed  in  section  8. 

Finally,  for  the  measurements  to  achieve  a  high  level  of  accuracy,  a  considerable 
amount  of  time  may  be  required  if  the  random  error  in  the  estimate  intensity  is  to 
be  small.  This  is  particularly  true  if  the  FFT  method  is  used. 
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13.  Conclusions 


A  system  has  been  developed  for  the  determination  of  the  sound  power  of  a 
source  based  on  the  principle  of  intensity  measurement.  In  this  system,  data  are 
taken  in  the  form  of  a  cross  power  spectnun  with  an  FFT  analyser.  The  data  are 
then  processed  using  a  micro-computer  for  intensity  calculation  and  synthesising 
of  third  octave  band  data.  By  analysing  sound  signals  in  the  frequency  domain^  a 
general  purpose  dual  channel  FFT  analyser  can  be  used  in  place  of  a  dedicated 
expensive  sound  intensity  analyser. 

This  report  has  demonstrated,  using  a  refereiKe  sound  source,  that  the  sound 
intensity  technique  can  be  applied  to  the  measurement  of  sound  power  with  a 
reasotuible  degree  of  accuracy  even  under  advert  conditioits  induding  the 
presence  of  background  noise.  Because  of  this  the  souitd  intensity  technique  can 
serve  as  a  useful  tool  to  identify  the  noise  emissions  of  RAN  surface  sh^  and 
subnnarines  in-situ. 
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15.  List  of  Symbols 


F  :  the  force  vector  in  space, 
a  :  the  acceleration  vector. 

Ej.  :  the  work  done  by  the  sound  field  in  the  r-direction. 

^r  :  distance  d  in  the  r-direction. 

A  :  area  of  the  surface  under  consideration. 

A  r  :  the  separation  between  2  points  in  space.  For  the  sound 
intensity  probe  it  is  the  separation  between  its  2 
microphones. 

t  :  the  amount  of  time  in  seconds. 

p  :  the  sound  pressure. 

Pr  ;  the  estimator  of  p  in  the  r-direction. 

PaiPb  •  pressures  measured  at  the  2  microphoites  of  the 

intensity  probe. 

Lp  :  the  sound  pressure  level  ref.  20  pPa 

Lp(90) :  the  sound  pressure  level  when  the  sound  wave  is  at  90“  to  the  probe 
axis. 

u  :  the  air  particle  verlodty  in  space. 

Uf  :  component  of  air  particle  verlodty  in  the  r-direction. 

Af  ;  the  estimator  of  Uj.. 

I  :  the  sound  intensity  vector  in  space. 

If  component  of  I  in  the  r-direction. 

:  the  estimator  of  1,.. 

L}  :  the  measured  intensity  level  ref.  1  pW. 

Ljres  :  the  residual  intensity  level  ref.  1  pW. 

4res(90)  •  residual  intensity  level  when  the  sound  wave  is  at  90“  to  the  probe 
axis. 
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Gab(  f )  •  the  cross  power  spectrum  between  the  2  microphones  of 

the  intensity  probe. 

Im  [  Ggb(  O  ]  :  imaginary  part  of  the  G^ib- 
p  ;  the  density  of  the  air. 

Af  :  frequency  increament  (resolution)  of  the  cross  power 
spectrum. 

N  :  Number  of  spectral  lines  in  the  cross  power  spectrum  G^b- 
k  :  the  wave  number. 

X  :  the  wave  length  of  the  sound  wave. 

<I>ab  :  the  detected  phase  difference  between  two  microphones  of  the  probe. 
5p|  :  the  pressure-intensity  index. 

SpiQ  ;  the  pressure-residual  intensity  index. 

L(]  ;  the  system  dynamic  capability. 


Intensity  p«>be, 


aerophone  spacer 


3520  with '12 


Intensity 


B.  The  calibration  of  the  phase  part  sumlied  with  the  microvhone  vair 


Serial  No.  1478121 
Resolution  0.05  dej>. 


Frequency 

(Hz) 


Phase  difference 


C.  The  sound  intensity  measurement  system: 


Figure  13:  The  sound  intensity  measurement  system.  Showing  from  left  to  right: 
the  CF-350  FFT  analyser,  the  intensity  probe  type  3520  and  the  microcomputer. 
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Appendix  2 

2 A.  Sound  source  twe  4205.  85  dB  nominated  outvut  - 125  Hz  band: 


2A.1  Full  octave  band  analysis  of  the  measured  output  power,  pressure-intensity 
index  and  pressure-residual  intensity  index: 


Microphone  spacer  At  =  12 mm 
Spectrum  fieqaency  nuige:  0  -  200  Hz 


Octave 

Band 


Freq 


Measttied  Power 


Mag  dB 

5.52E-05  7%B2 

5.16E-04  87.13 


2B.  Sound  source  tuve  4205.  85  dB  nominated  output  - 125  Hz  baud: 


2B.1  Full  octave  band  analysis  of  the  measured  output  power,  prcssure-intensit 
index  and  pressure-residual  intensity  index: 


Miciophone  spacer  A  r 
Spectrum  fnquency  lange 


Measuied  Power 


dB 

5.25E-05  77.2 

4.27E-04  86.304 

1.12E4B  703 


Below  are  the  results  as  displayed  by  the  sound  intensity  program. 
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2C.1  Full  octave  band  analYwis  of  the  measured  output  power.  Dressure-intgnsii 


index  and  pressure-residual  intensity  index: 


Microphone  spacer  A  r  s  12  mm 
Spectrum  frequency  lange:  0-500  Hz 


Octave 

Band 

Measured  Power 

Spl 

5pI0 

Freq 

Mag 

dB 

dB 

dB 

31.5 

2.50E-07 

53.981 

7.147 

7.253 

63 

1.57E-06 

61.966 

6.081 

9.516 

125 

6.45E-05 

78.098 

2.899 

13.426 

250 

3.20E-04 

85.0577 

2.939 

17.526 

500 

2.25E-05 

73.513 

2.926 

16.932 

2C.2  Third  octave  band  analysis  of  the  measured  output  power 

Below  are  the  results  as  displayed  by  the  sound  intensity  program. 
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2D.  Sound  source  type  4205.  85  dB  nominated  output  -  500  Hz  band: 

2D.1  Full  octave  band  analysis  of  the  measured  output  power,  pressttre-intensit 
index  jnd  pressure-residual  intensity  index: 


MJerophone  spacer  A  r  =  12  mm 
Spectrum  frequency  range:  0-1  kHz 


Octave 

Band 

Measured  Power 

8pl 

®pro 

Freq 

Mag 

dB 

dB 

315 

6.51E-07 

65.63 

5.857 

6.9 

63 

2.32E-07 

53.659 

7.391 

9.765 

125 

S.97EW 

59528 

5.264 

12.558 

250 

1.62E-05 

72.082 

2.882 

16.954 

500 

2.89E-04 

84.615 

3.248 

18.171 

l.OOE+03 

5.94E^ 

77.735 

3.542 

22.03 

2D.2  Third  octave  band  analysis  of  the  measured  output  power: 

Below  are  the  results  as  displayed  by  the  sound  intensity  program. 
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gllK^  JjgTiiEffI  n 


2E.1  F 

index  and  Diessure-residual  intensity  index; 


T*  ■>Tn  niTH  ^  > 


ressure-intensi 


Microphone  spacer  A  r  *  12  mm 
Spectrum ftc4|aency range:  0-2kHz 


Octave 

Band 

Measured  Power 

»pl 

®pI0 

Frecj 

Mag 

dB 

dB 

dB 

63 

2.09E-06 

58.378 

8.159 

8.839 

125 

1.31E-07 

51.17 

9.424 

11.978 

250 

2.00E-07 

53.02 

6.955 

16.027 

500 

2.72E-05 

74.339 

3.115 

19.121 

l.OOE+03 

2.90E-04 

84.62 

3.152 

20.084 

2.00E+03 

1.07E-(B 

70.313 

2.471 

20.888 

Below  are  the  results  as  displayed  by  the  sound  intensity  program. 


2F.  Sound  source 


type  4205.  85  dB  nominated  output  -  2K  Hz  band: 

2F.1  Full  CKtave  band  analysis  of  the  measured  output  power,  pressure-intensity 
index  and  pressure-residual  intensity  index: 


Mkfophone  spacer  A  r  -  12  mm 
Spectrum  frequency  range:  0-5KHz 


Octave 

Band 

Measured  Power 

8pl 

8pI0 

Fretf 

hAag 

dB 

dB 

dB 

250 

3.25E-08 

45.113 

11.518 

14.869 

500 

3.33E-08 

45.23 

9.324 

18.342 

l.OOE+03 

3.56E-05 

75.511 

3.084 

20.65 

2.00E+03 

2.%E-04 

84.72 

2.988 

22.733 

4.00E+03 

2.75E-05 

74.399 

3.472 

21.445 

2F.2  Third  octave  band  analysis  of  the  measured  output  power 

Below  are  the  results  as  displayed  by  the  soiuid  intensity  program. 
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2G.  Sound  source  type  4205.  85  dB  nominated  output  -  4K  Hz  baud: 

2G.1  Full  octave  band  analysis  of  the  measured  output  power,  pressttre-intensitv 
index  and  pressure-residual  intensity  index: 


Microphone  spacer  A  r  =  12  mm 
Spectrum  frequency  nuige:  0  •  lOK  Hz 


Octave 

Band 

Measnnd  Power 

8pl 

S 

IC 

Freq 

Mag 

dB 

dB 

dB 

500 

2.29E-09 

33.597 

10.654 

18.657 

l.OOE+03 

5.75E-07 

57.6 

3.774 

21.219 

2.00E+03 

2.10E-05 

73.36 

3.031 

22.375 

4.00E+03 

2.57E-04 

84.14 

3.145 

20.019 

8.00E+03 

l.%E-05 

72.9 

4.022 

17.827 

2G.2  Third  octave  band  analysis  of  the  measured  output  power: 

Below  are  the  results  as  displayed  by  the  sound  intensity  program. 
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